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ABSTRACT 


Consider  fractional  factorial  designs  of  resolution  IV,  2., a.  where  we  wish 
to  estimate  only  ine  main  effects  but  the  2-  factor  interactions  are  not 
negligible.  Such  designs  with  desirable  size  ere  greatly  needed  both  in 
agriculture  and  industry,  and  both  in  uniresponse  and  multiresponse  experiments. 
The  usual  completely  orthogonal  designs  irr’olve  N  runs,  where  I!  is  a  ..-uitiple  of 
8.  In  many  situations,  we  have  a  set  of  exactly  N  homogeneous  experiment  units, 
where  N  is  not  divisible  by  S,  For  example,  we  may  have  H  =  22  new  jet  bombers 
s*f  •'  certain  kind  being  developed  for  defense  purposes.  Here  the  sampling  units, 
namely  the  ')ets,  are  expensive,  and  furthermore  cannot  easily  be  increased  in 
number  ju.it  for  the  sake  of  our  experiment.  If  we  want  to  test  these  jets  under 
a  factorial  design  of  resolution  IV  using  the  present  available  designs,  then  we 
can  use  oalv  10  ot  them,  since  16  is  the  multiple  of  8  nearest  to  but  not  greater 
than  22.  This  voxsjJ  result  in  a  loss  of  6/22  of  the  available  information. 

Tnus  the  purp ..  e  of  this  paper  is  to  obtain  good  nonortnogonal  or  irregular 

•  HI  • 

designs  of  resolution  7X  for  the  2  series.  Besides  being  of  desirable  size,  a 
nonorthogonal  design  should  be  good  with  inspect  to  its  covariance  matrix  V  of 
the  estimates.  In  this  paper,  such  designs  (with  even  N)  are  obtained.  These 
designs  are  optimal  with  respect  to  the  trace ,  determinant  and  the  largest  root 
criteria  which  are  shown  to  b:  equi'. uicnt.  In  other  words,  amor.g  all  possible 
designs  a  given  value  of  N,  our  designs  minimize  the  trace,  the  determinant  and 
the  largest  root  of  V. 
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Summary .  In  this  paper,  we  develop  a  general  theory  of  balanced  2  fractional 

.  1  -  '  r  4. 

factorial  designs  which  permit  estimation  of  main  effects  orthogonally  to  2- 
factor  interactions  and  the  general  mean,  whose  size  N  is  desirably  small,  and 
which  are  optimal  with  respect  to  various  standard  criteria  involving  the  variance- 
covariance  matrix  of  the  estimates.  For  various  practical  values  of  m  and  N, 
a  method  is  given  by  which  such  optimal  designs  can  be  easily  obtained  from  known 
balanced  incomplete  block  designs,  (BIBD's). 

1.  Introduction 

Fractional  factorial  designs,  discovered  by  Finney  (1945),  are  finding 
increasing  use  in  biological,  industrial,  and  defense  research.  A  large  number 
of  authors  have  contributed  to  the  development  of  the  various  facets  of  the  theory. 
For  the  reader  interested  in  an  introduction  to  the  same ,  a  list  of  such  authors 
is  included  in  the  references  at  the  end.  We  stress  that  this  list  is  by  no  means 
exhaustive;  it  is  only  illustrative. 

In  many  cases,  the  main  effects  are  of  immediate  or  primary  interest,  but 
ti.e  two-factor  interactions  cannot  be  assumed  negligible.  "Good"  designs  for 
various  symmetrical  and  asymmetrical  factorials  need  therefore  to  be  developed. 

The  properties  of  "goodness"  that  are  most  important ,  and  most  often  even  necessary 
include  the  following:  (a)  The  design  should  be  of  small  size,  i.e.,  the  number 
of  observations  should  be  the  desirable  minimum,  (b)  The  design  should  have 
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balance  or  at  least  partial  balance.  This  lends  facility  not  only  to  the  analysis 
of  the  design,  but  also  to  the  interpretation  and  understanding  of  the  results. 

(c)  The  design  should  further  satisfy  some  reasonable  optimality  condition  on 
the  variance  covariance-matrix  of  the  estimates. 

Let  L  denote  tne  (best  linear  unbiased)  estimate  of  the  parameter  vector  L^. 

The  usual  variance  criteria  under  (c)  above  include  (i)  trace,  (ii)  determinant, 
and  (iii)  largest  root  criterion.  A  design  T  is  said  to  satisfy  these  criteria  if, 
respective Ly,  the  trace,  determinant,  or  largest  root  of  is  a  minimum  for  this 
design  within  the  class  of  all  possible  designs  of  fixed  size.  Here  =  Cov  (Z) , 
when  L  is  estimated  using  T. 

In  this  paper  we  shall  discuss  designs  satisfying  all  the  criteria  mentioned 
above. 

2 .  Some  Preliminaries 

We  shall  use  a  notation  similar  to  the  one  in  Bose  and  Srivastava  (19641, 
and  for  the  facility  of  the  reader,  repeat  (without  detailed  proof)  some  basic 
theory  of  fractional  designs  developed  therein.  Treatment  combinations  (and  also 

h  ^m 

their  "true  effects")  are  denoted  by  (a^  a 2  . ..,  a  *}  ,  or  equivalently 

(j.,  •  ••>  j  )  where  i's  e  (0,1).  The  k  factor  interactions  between  factors 

j.  z  m 

A.  ,  A.  ,  ...,  A.  ,  iay  is  denoted  by  (A.  A.  ,  ....  A.  ).  As  usual,  the  interactions 
11  ”2  1k  X1  12  '  xk 

are  defined,  in  symbolic  notation  by 


(2.1)  AxA2  ...  Ak 


,  1  0.,  1  0.  ,  1  0.,  1  0  .  ,1  0, 

(a1-a1)(a2~a2)  ...  <ak-ak)(ak+1+ak+1) . . .  (am  taj. 


and  similarly  for  any  set  of  k  factors,  where  1  <  k  <  n.  The  symbol  <f>,  which 


denotes  the  grand  total  of  all  2'  assemblies,  is  also  given  by  (2.1),  if  the 


r.h.s.  is  assumed  to  have  plus  sign  in  each  bracket. 
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Let  £  denote  the  column  vector  of  A’s  in  the  standard  order: 

(2.2)  £*  =  [♦;  Am;  A^,  A^,...,  A^Ag,  ...] 

If  £  (2m  x  1)  denotes  the  vector  of  assemblies,  then  (2.x)  implies  £  =  D_r_, 

where  D  is  a  (2m  x  2m)  matrix,  and  the  sum  of  products  of  the  corresponding 

-m/2 

elements  in  any  two  rows  of  D  is  zero.  It  is  easily  seen  that  2  D  is  an 
orthogonal  matrix  and 

(2.3)  £  =  2-InD'£  . 

Let  £  be  partitioned  as 

(2.4)  £'  =  (£' :I^) 

where  £  (say,  (v  x  l),  where  v  -  1  +  m(m+l)/2)  is  the  column  vector  containing 

all  effects  up  to  and  including  all  the  two  factor  interactions,  and  I^(with 

(2m  -  v)  elements)  is  the  vector  containing  all  3- factor  and  higher  order  interactions. 

In  this  paper,  is  assumed  zero.  Then  (2.4)  implies 


(2.5) 


t  =  2  mD!L 
—  u — 


\  •  *1 
where  is  the  matrix  obtained  by  cutting  out  the  (2  -  v)  columns  of  D’  which 

correspond  to  1^  in  (2.3). 

Let  T  be  a  fractional  design,  i.e.  a  set  of  assemblies  in  which  any  given 
treatment  may  not  occur  or  may  occur  once  or  more  times.  Let  the  expected  values 
of  the  assemblies  in  T,  written  in  the  form  of  a  vector,  be  denoted  by  z*  where, 
it  is  assumed  there  are  no  block  effects.  Let  E'  be  the  matrix 
obtained  from  by  cutting  out  the  last  (2m  -  v)  columns  corresponding  to  1^, 
and  also  by  omitting  (or  repeating)  the  rows  corresponding  to  treatment  combinations 
omitced  (or  repeated)  from  t.  to  get  z*.  Note  that  the  rows  of  E'  are  arranged  in 
such  a  way  as  co  correspond  to  the  elements  of  £*.  Let  z  be  the  vector  of 
observations  corresponding  to  z*.  Assuming  no  block  effects  the  model  is 


4 


(2.6)  Exp  (z)  =  z*  =  2"mE'L 

var  (z)  =  c2i 
—  n 

where  °2  is  unknown,  and  I  is  the  (n  *  n)  identity  matrix. 

Then  the  normal  equations  are 

(2.7)  EE'£  =  Ez 

•  where  o  stands  for  the  estimate  of  =  2  L^. 

Let  T  be  a  fractional  design  containing  N  assemblies.  Then  the  symbol 


u1/  r 
il12*  *  *  Xr 


W  3l  32  \ 

or  X(a.  a.  ...  a.  ) 

It  lo  1 

12  r 


denotes  the  number  of  assemblies  as  such  that  the  symbol 


D  t  3  2  j  r  *i 

a.  a.  ...  a.  occurs  as  a  part  of  os.  The  symbols  a.  are  to  be  regarded  as 

X1  12  xv  1 


indeterminates  and,  as  such,  expressions  like  the  above  are  operated  upon,  over 
the  field  of  real  numbers,  like  ordinary  products  of  indeterminates. 

Consider  the  set  P  of  all  polynomials  with  real  coefficients  in  the  symbols 
a^  ,  (1  <  i  <m,  j  =  C,  1)  and  of  degree  1  in  each  symbol,  such  that  no  term 
involves  both  a?  and  a^,  for  any  i.  In  the  sequel,  we  shall  be  ;oncerned  with 
only  such  polynomials;  more  information  on  these  can  be  had  from  Srivastava  (1967). 
Then  we  define 


(2.8)  UBp1  +  6p2)  =  8X(p1)  +  5X(p2), 

for  all  p1 ,  p2  c  P.  Thus  for  example  A(2a'}a2  -  3a^a^)  =  2A^2  -  3A^. 

If  T,  and  T2  are  two  fractional  designs,  then  T^  @  T2  shall  denote  the  design 
in  which  any  treatment  combination  u>  occurs  (r^  +  r2)  times,  provided  u>  occurs 
respectively  r^  and  r2  times  in  T^  and  T2-  Then  if  r  P,  we  have 

(2.9)  3.(p,T^$T2 )  =  X(p,T)  +  A(p,T2), 

where  A(p,T)  denotes  the  value  of  A(p)  obtained  for  the  set  of  assemblies  T. 


>  *va»rerar.;« htcsu.wv-.ii^wIw  r^v<r.-»TO-»rwnr--  ^pro'WwnwiwwraM 
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Given  any  fraction  T  the  matrix  EE'  can  be  directly  expressed  using  the 

X-operator.  Every  row,  and  hence  every  column,  of  EE'  corresponds  to  exactly 

one  element  of  L^.  Indeed,  the  element  in  the  (i,j)  position  of  EL'  corresponds 

th  vh 

to  the  elements  in  the  i  and  i'-1  position  of  L.  An  element  of  can  be  expressed 
jl  j2 

as  A.  A.  where  j  ,  j  =  0  cr  1.  Then  corresDonding  to  {p}  (the  general  mean), 

11  x2 

(A.)  (the  main  effects),  and  {A. A.)  (the  set  of  two-factor  interactions,  the 
i  i  j 

matrix  (EE1 )  can  be  nartitioned  in  the  form 


(p)  (A.)  {A. A.} 

i  l  3 


(2.10) 


EE’  = 


M11  H12 


22  23 


(y) 

<  <v 

M,_  (A. A.: 

33  J  13 


Theorem  2.1.  We  have,  for  distinct  i,  j,  k,  and  £, 

(2.11)  e(p)  =  e(A.,A.)  =  e(A. A . ,  A. A.)  =  N 

11  13’  13 

(2.12)  c(A.,A.)  =  e(A.  ,A. A . )  =  X ( a"!"  -  a?) 

13  113  3  3 

(2.13)  e(A.,A.)  =  c(p,A.A. )  =  e(A.A  ,  A. A.)  =  X(a*  -  a?)(aA  -  a?) 

13  13-  1  k’  3  k  1  r  3  3 

(2.14)  e( A.,A.A  )  =  X(a'l'  -  a?)(a‘!'  -  a?)(a^  -  a®) 

13k  1  13  3  k  k 

(2.15)  e(A.A.,A  A  )  =  X(a'l‘  -  a?)(at  -  a?)(a^  -  a°)(a„  -  a°) 

.t  3  k  £  1  13  3k  k  £  £ 

The  value  of  the  X-operator  for  the  above  polynomials  can  be  calculated  by 
using  the  definition  at  (2.6).  For  example 

(2.16)  e(A.  ,A.A,  )  =  \(a}  -  a?)(al  -  a?)(a?;  -  a^) 

13k  1  1  3  3  k  k 

,,111  110  101  011  100  010 

=  XCa.a.a,  -  a. a. a,  -  a. a. a,  -  a. a. a,  +  a. a. a,  +  a. a. a, 

1  3  y  1  3  k  13k  1  j  k  1  3  k  1  3  x 

0  0  1  0  0  0 

t  a. a. a,  —  a. a. a, 

13k  13k 

_  ,111  ,110  ,101  ,011  ,  :oo  ,010  ,ooi  ^  ,000 

ijk  ijk  ijk  13  k  13k  13k  13k  ijk 
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3,  The  2"  Factorial  Experiment 

Let  T  be  a  fractional  design  for  a  2m  factorial  experiment  with  1  assemblies. 
Recall  that  T  is  regarded  as  an  (n  x  N)  matrix  of  0's  and  i'r  where  each  column 
corresponds  to  an  assembly. 

Pef.  3.1.  (a)  T  is  said  to  be  of  resolution  IV,  if  ail  main  affects  are  estimable 

under  the  assumption  that  3-factor  and  higher  order  interactions  are  zero.  (Notice 
that  we  have  no  such  assumption  regarding  2- factor  interactions ;  they  might  possibly 
be  nonzero).  If,  further-,  the  estimates  of  the  main  effects  have  zero  correlation 
with  the  estimates  of  every  estimable  linear-  combination  of  the  general  lean  y  anc 
the  two- factor  Interactions,  then  T  is  said  to  be  cf  resolution  JV*. 

(b)  In  general,  in  any  given  situation,  if  T  is  such  that  ail  parameters 

of  interest  are  estimable,  then  T  is  called  nonsingular.  Thus,  a*  this  paper,  the 
properties  of  nonsingularity  and  resolution  IV  are  equivalent. 

(c)  The  matrix  T  (re  x  K)  is  said  to  be  (1,0)  symmetric  of  strength  t  (j^m), 
if  every  (t  x  N)  subroatrix  remains  unchanged  by  interchanging  the  symbols  0  and  1, 
except  for  the  order  in  vhich  the  columns  appear.  If  t  -  m,  T  is  called 
"completely  (1,0)  syrisnetric".  (Clearly,  strength  t  implies  strength  i-,  for 

r  <_  t.  Also,  the  authors  havs  examples  of  arrays  T,  which  are  of  e  certain 
strength  r,  but  not  of  strength  t,  for  some  t  >  j). 

(d)  T  is  said  to  be  balanced  if  ccv  (y,  A.),  var  (A.),  and  cov  (A.,  A.) 

1  i  s  j 

(i  f  j),  are  independent  of  i  and  j. 

Theorem  3.1.  T  is  of  resolution  IV*  if  and  only  if  it  is  (1,0)  syumetric  of 
strength  3,  and  (EE1)^  is  nonsingular. 

Proof.  That  (EE')  be  nonsingular  is  obviously  necessary.  Ma5.n  erfeccs  will  be 
orthogonal  to  general  mean  and  two  factor  interactions  if  and  only  if  and 
in  equation  (2.10)  are  zero  matrices.  Thus  we  have  from  ( 2. 11 ) - ( 2.15) : 

(3.1)  0  =  e(u,A.)  =  A^  -  A?  =  e(A.,A.A.);  j  =  1,  2,  m. 

]  ]  3  113  J  ’’  ’ 
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Proof . 

Suppose  X^  =  (u ,  for  all  i  /  j. 

Now 

(3.6) 

e(A.  ,A. )  =  X(al‘ 
11  1 

+  a?)  =  N; 

1 

(3.7) 

e(A. ,A . )  =  X(a? 
13  1 

0W  1  0, 

-  a. )(a.  -  a. ) 
i-  3  3 

11  10  01  00  vr 

-  X . .  —  X -  .  —  X..  +  X..  =  —  N > 

13  13  13  13 

Hence 

(3.8) 

=  (2N  -  4u)I  +  (4oj  -  N)J 

22  y  m  mm 

;  and 

(3.9) 

»22  *  <8  -  h)i„ 

+  hJ  ,  where 

fissn 

(3.10a) 

1 

m-1 

®  m[4(m-l)o)  - 

(m-2jKT  " 

m(2N-4oj) 

(3.10b) 

U  -  A 

1 

m[4(m-l)(o  - 

(m-2)N] 

m(2N  -  4a.) 

(3.10c) 

g  -  h  =  1/C2N  - 

4u] 

Hence, 

if  xH  is  a  constant  u>  for  all  i 
ij 

i  j,  then  T  is  balanced.  On  the  other  hand 

if  T  is  balanced,  ther.  has  all  diagonal  elements  same  and  all  off-diagonal 
element  same,  so  that  is  of  the  form  (3.8)  for  some  value  of  u>.  Hence  (3.6) 
holds  for  all  i  /  j.  This  completes  the  proof. 


It  follows  from  the  above  theorem  that  if  T  is  a  balanced  fractional  design, 

tnen  T  is  the  incidence  matrix  of  a  BIBDU,  (i.e.  a  BIBD  with  possibly  unequal 

block  sizes),  and  h'  ing  parameters  of  the  form  (v*  =  m,  b*  =  N,  r*  =  N/2,  and 

X*  =  uj).  This  observation  provides  a  simple  method  for  obtaining  balanced  fractiona 

designs.  For  example,  if  T^  is  the  incidence  matrix  of  a  BIBDU  with  parar ;ters 

(v*  =  m,  b  =  N/2,  r*,  X*)5  then  T  =  T^  ®  is  a  balanced  fractional  design  where 

T  is  the  (1,0)  complement  of  T^ ,  i.e.  T^  is  obtained  from  ^  by  interchanging 

.  .  th 

1  and  0.  It  is  obvious  that  T  is  completely  (1,0)  symmetric.  Consider  the  i 

and  jth  rows  of  T^.  Then,  for  T..  ,  xH  =  X*  and  X??  =  b*  -  2r*  +  X*,  and  by  the 

method  of  construction  of  T,  it  is  clear  that  X^(T)  =  xH(T  }  +  X?.(T1).  Hence, 

13  13  1  1 


for  1  we  have 
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(3.11)  0)  =  b*  -  2(r*  -  X) 

A  BIBDU  may  be  easily  obtained  from  an  ordinary  BIBD,  by  cutting  out  some  of  the 
treatments  and  perhaps  adjoining  some  blocks  containing  all  treatments  and/or 

/I  »t 

some  blocks  containing  no  treatments.  As  an  example  take  a  (BIBD),  with  parameters 
(v*,  b*,  r*,  k*,  X*)  and  cut  out  <  v*  varieties.  Next  adjoin  b^  blocks  containing 
each  of  the  remaining  (v*-  v  )  varieties,  and  b^  blocks  containing  no  varieties. 

The  resulting  EIBDU  has  parameters  ((v*  -  v^) ,  (b*  +  b^  +  b0),  (r*  +  b, ) , 

(X*  +  b  )).  (It  may  be  beneficial  to  the  interested  reader  to  remark  here  that 
the  incidence  matrix  of  a  BIBDU  is  a  partially  balanced  (PE)  array  of  strength 
2,  the  latter  being. defined  for  example  in  Chakrfivarti  (1963)  or  Srivastava  (1967  )). 

We  now  consider  the  non-singularity  and-  optimality  of  the  balanced  fractional 

i 

designs  obtained  in  this  manner. 

\ 

x  l 

Lemma  3.1.  Let  T  be  a  balanced  design  as  in  Theorem  3.2.  Thai (a)  The  characteristic 

roots  of  M22  are  [4(m-l)w  -  (m-2)N]  and  [2N  -  4w],  with  multiplicities  1  and  (m-1) 

1 

respectively, 

(b)  Tfe  characteristic  roots  of  M22  are  [g  +  (m-l)h]  and  (g-h)  with  multiplicities 
1  and  (m-1)  respectively.  Also  tr(M22)  -  mg  and  | I  =  Eg  +  (m-l)'n](g-h)m  1  where 
g  and  h  are  given  by  (3.10). 

Proof:  Let  k,  and  k„  be  constants.  Then  the  matrix  k.I  has  a  root  k,  wich 
12  1ml 

multiplicity  m,  ard  the  matrix  k2Jmm  has  roots  mk2  and  0  with  multiplicity  1  and 
(m-1)  respectively.  Since  I  and  J  consnute,  the  roots  of  (k,I  +  k„J  )  are 
(k^  +  mk2)  and  k^  with  multiplicities  1  and  (m-1).  The  proof  is  completed  by 
appropriate  substitution  for  k^  and  k2,  a^d  noting  that  the  trace  and  determinant 
are  respectively  the  sum  and  the  product  of  the  roots. 

Theorem  3.3.  A  necessary  and  sufficient  condition  that  a  balanced  fractional 


design  T  be  nonsingular  is  that 


11 

(3.12)  C(m-2)/4(m-l)]|T|  <  <o(T)  <  ( 1/2) j i j  , 

where  T  denotes  the  number  of  assemblies  in  T,  and  m(T)  is  the  value  of  a)  associated 
with  T. 

Proof:  The  matrix  M  =  EE*  is  obviously  positive  semidefinite .  Thus  X 

*./.  ±  i  22 

is  nonsingular  if  and  only  if  ail  roots  of  MQ2  are  positive  and  tne  result  folxows 
from  lemma  3.1. 

Corollary  3.1.  Let  m  be  even,  and  let  T  be  the  incidence  matrix  of  a  BIBD  with 
k*  =  m/2,  ar.d  v*  =  m.  Then  T,  considered  as  a  fractional  design  is,  singular. 
However,  the  balanced  design  obtained  by  adjoining  to  T  (at  least)  one  block 
containing  all  varieties  and  (at  least)  one  Dlock  containing  no  variety,  is 
nonsingular. 

Proof.  Since  v*r*  =  b*k*,  and  r*(k*-l)  =  X*(v*-1),  we  get  under  the  stated 
conditions:  r*  =  b*/2,  |?|  =  b*,  and  X*  =  w.  Thus  we  have  uj  =  b*(m-2)/4(m-.l) . 

Hence,  by  (3.12),  T  is  singular.  The  fractional  design  T*  obtained  by  adjoining 
the  two  blocks  has  | T* |  =  b*  +  2,  and  o>(T*)  =  u  +  1.  The  proof  is  easily 
completed  by  checking  that  these  values  satisfy  (3.12). 

Corollary  3.2.  Let  T  =  T,  @  1^ ,  where  is  the  incidence  matrix  of  a  BIBOU 
with  parameters  v*  =  m,b*,r*,  X*.  Then  T  is  nonsingular  if  and  only  if 
(3.13)  0  <  r*  -  X*  <  mb*/4(m-l). 

Proof.  Clearly  |t|  =  2b*,  and  by  (3.11),  u  =  b*  -  2(r*  -  X*).  Then  (3.12)  gives 
2b*(m-2)/4(m-l)  <  b*  -  2(r*  -  X*)  <  2b/2 
which  leads  tc  (3.13). 

Let  denote  the  class  of  all  balanced  nonsingular  completely  (1,0)  symmetric 
fractional  designs.  Then  N  is  even,  since  X^  =  X°  for  all  T  Let/^  denote 

the  class  of  all  T  z'f,  with  |t|  -  N. 
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Theorem  3,4.  Let  T  s/f,.  Then  N  _>  2m,  and  ng  (the  number  of  degrees  of  freedom 

for  error)  satisfies 

(3.14)  (N-2m)/2  <  n  <  N-2m. 

—  e  — 

Proof.  Let  c  =  4m  -  N.  Then,  M.„  =  (N-c)I  +  cJ  .  Let  M*  be  the  principal 
-  22  m  mm 

submatrix  of  (EE')  having  the  rows  and  columns  corresponding  to  {p  ,A  A  ,A  A0 , . . .  ,A  A  }. 
Then  it  can  be  checked  by  direct  calculation  that  =  M*.  Also,  since  T 
T  and  hence  is  nonsingular.  Since  and  are  obviously  zero  matrices, 
this  implies  Rank  (EE-')  >  2m.  But  Rank  (EE')  =  Rank  (E)  _<  N;  hence  N  >_  2m. 

Also,  the  number  of  degrees  of  freedom  for  error  is  ng  -  N-Rank  (M)  and  thus 
n  <_  N  -  2m.  finally,  let  T  e  and  let  (M**)  be  the  principal  submatrix  of 
(M)_  corresponding  to  p,  and  2--factor  interactions.  Now  T  can  be  divided  into 
two  pa-V  7^  and  T 2  =  T^.  Then  it  can  be  checked  that  M**  is  the  came  for  T^ 
as  for  ,  except  for  a  constant  multiplier.  Thus  R(M**)T  =  R(M**)T  <_  N/2,  and 

R(M)T  =  R(^22)t  +  I  m  +  N/2 

Hence  n  >  N  -  (N  +  2m)/2  =  (N  -  2m)/2,  proving  (3.14). 

Theorem  3,3.  net  T  e^.,  then  T  is  trace  optimal  in ^  if  one  of  the  following 
holds  for  sore  positive  integer  a. 

'3.17)  ( i )  N  =  4c  and  u  =  u 

(3.18)  (iij  ;<  -  4ct  +  2  and  iu  =  a  +  1 

^roof :  From  lemrr:c.  3.1,  and  (3.8)  we  have 


trtM"^)  =  *-g 


m-1 


—  + 


[4(m-l)«  -  (m-2-)Nj 


f(«). 


2N-4m 

say,  where  f(m)  is  a  ft  notion  of  m  for  fixed  m  and  N.  Suppose  “  were  a  continuous 
variable.  Then 

df(aj) 


(3.19)  f'(to)  = 


du 


-4m(m-l)(N-4u)[4h)(m-2)  -  (m-4)N] 
[2N-4u]^[4(m-l)u)  -  (m-2)N]-'- 
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and  f (w)  =  0,  when  w  =  N/4  or  to  =  (m-4)N/4(rn-2) .  The  root  cu  =  (m-4)N/4)(ra~2) 
is  extraneous  since  by  Theorem  3.3.  us  >  N(m-2)/4(m-l) ,  Since  the  denominator  of 
(3.19)  is  always  positive  it  is  easy  to  see  that  f * (to)  <  0  if  w  <  N/ 4  and  >  0 
if  w  >  N/4.  Hence  f(w)  attains  an  absolute  minimum  when  u)  =  N/4  and  (3.17) 
is  proven.  Now,  suppose  N  =  4a  +  2.  Since  to  can  take  only  integer  values, 
to  =  N/4  is  not  possible.  Hence  we  must  find  the  integer  value  of  to  which  makes 
f(to)  as  small  as  possible.  Since  f(to)  is  strictly  decreasing  on  (N(m-2)/4(ir.-l)  ,N/4) 
and  increasing  on  (N/4,N/2)  the  absolute  minimum  can  occur  either  at  to  =a  =  N/4  -  1/2 
or  to  =  a  +  1  =  N/4  +  1/2.  But  it  can  be  checked  that 


e [  \  _  16m(m-l)( 2-m)  .  „ 

na  +  i.)  -  net;  -  (N_2)(N4?)rN+2(n!-l)][N-2(m-l)3  -  ’ 


since  m  >_  2,  and  N  >_  2m.  This  completes  the  proof. 

Theorem  3.6.  Let  T  e  then  if  T  is  trace  optimal  in  2^,  if  T  is  determinant 

optimal  and  conversely. 

Proof.  Let 


f^Coi)  =  | Mp 2 1  =  [N(m-l)<o  -  (m-2)N][2N  -  4a)]"1  1, 

We  must  show  that  for  any  (even)  N.  Im^jI  and  tHM^)  become  minimum  for  the  same 
value  of  (o. 

Then 


f|(<o)  =  4m(m-l/(2N-4u))m  ^(N-4u>), 

The  roots  of  fj^u>)  =  0  are  oj  =  N/2  or  e>  =  N/4.  Among  these  m  =  V/2  is  not  ir  the 
domain  of  f^.  Again  ic  is  easy  to  see  that  f^(o>)  is  increasing  on  (N(m-2)/4(m-l) « 
N/4)  and  decreasing  on  (N/4,  N/2);  thus  an  absolute  maximum  is  reached  at  a  =  n/4 
and  the  theorem  is  proven  if  N  =  4a  ,  a  an  integer.  When  N  =  4a  +  2,  we  have 
f  (a  +  1)  -  f.  (a)  =  (N-2)m  -  (N+2)m  +  2m[(N-2)I!,"i+  (N+2)m"1] 

1  j. 

=  (N+2)m  1Q  >  say,  where 
m'  J 

Qin  =  (N+2)m-1[2m( l+S1"-^)  -  4(1-  g)_1(l-Bin)3, 
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where  3  =  (N-2)/(N+2).  Now  =  32(N+2)  >  0.  This  proves  the  result  for  m=3, 

since  f^ot  +  1)  -  f^a;  >  0.  For  larger  n.,  it  is  clearly  sufficient  to  show 
that  0  is  an  increasing  function  of  m.  But  it  can  be  checked  that 

-  Qm_1  =  2(l-8)[(l-8m"2)  +  6(l-8rn_3)  +  ....  +  6m_2(l-8m'ni)], 
which  is  ^  0,  since  8  <  1.  This  completes  the  proof. 

Theorem  3.7.  Let  T  then  trace  optimality  of  T  implies  maximum  root  optimality, 

and  conversely. 

Proof:  From  lemma  3.1  the  characteristic  roots  of  ate  [4(m-l)w  -  (m-2)N]_1 
and  (2N~4to)  \  Now 

(3.20)  [4(m-l)to  -  (m-2)N]  < ,  =,  or  >  2N-4  ,  according  as  to  <,  =,  >  N/4. 

If  N=4a,  it  is  obvious  from  (3.20)  that  the  maximum  root  of  takes  the  smallest 
value  when  to  =  N/4.  When  N  =  4a  +2,  the  max.  root  of  is  [4(m-l)to  -  (m-2)N]_1 
or  [2N  -  4ui]  1  according  as  to  =  a  or  a  +  1. 

But  (3.20)  shows  that  among  these,  the  value  of  the  max.  root  of  is  smaller 
at  u  =  a  4-  1.  This  completes  the  proof. 

Theorem  3.8.  Let  T  be  the  incidence  matrix  of  a  BIBDU  with  parameters 
(v*  =  m,b*,r*,A*)  such  that  T  =  T^  @  Tn  e  Let//  denote  the  class  of  fractions 

obtainable  from  BIBD’s  in  this  tanner.  Then  T  is  trace  optimal  within  Zt  if  there 

exists  an  integer  a  such  that  any  of  the  following  conditions  hold: 

(i)  b  =  4a  +0;  0=0,  1,  or  2,  and  r*  -  A*  =  a; 

(ii)  b  =  4a  +  3,  and  r*  -  >.*  =  a  +  1. 

Proof:  It  has  already  been  shown  that  if  T  =  T^  ©  T^,  then  to  =  b  -  2(r*  -  A*) 
and  N  =  2b*.  Thus  r*  -  A*  =  (b-to)/2. 

When  b*  =  4a,  we  gee  N  =  4(2a),  and  the  optimum  value  of  to  is  2a  so  that  r*  -  a*  =  a. 
When  b*  =  4a  +  1,  we  get  N  =  4(2a)  +  2,  and  the  optimum  value  of  to  is  (2a  +1), 
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so  that  r*  -  X*  =  a.  Next, 
algebraically  optimum  value 
which  is  not  possible  since 
possible  optimum  value  of  w 
proof  of  theorem  3.5.  Then 


suppose  b*  =  4a  +  2 .  Then  N  =  4(2a  +  1)  and  the 
of  a)  is  (2a  +  1)  which  implies  (r*  -  A”)  =  (a  +  1/2), 
r*.  A*  and  a  are  integers.  Hence  the  conbinatorially 
is  either  2a  or  (2a  +  2).  Define  f(m)  be  ss  in  the 


m 


f(2a  +  2)  -  f(2a)  =  J 


-64(m-l)(m-2) 


(|  -  2)[N+4(m-1 ' J(N/2  +  2)[N-4(m-l)] J 


which  is  <,  or  >  0,  acceding  as  N  >,  or  <  4(m-l).  But  N  <  4(m-l)  implies 
a)  <  N(m-2)/4(m-l).  Hence  w  -  2a  +  2.  Tnus  part  (i)  is  proven,  since 
r*  -  A*  =  [(4a  +  2)  -  ;2a  +  2)]/2  =  u. 

If  b*  =  4a  +  3,  then  N  4(2a  +  1)  +  2  and  from  theorem  3.5,  che  algebraically 
optimum  value  for  «j  is  (2a  +  2).  However,  w  =  2a  +  2  leads  to  r*  -  A*  -  (a  -r  1/2). 
Thus  the  choice  for  actual  optimum  is  between  u  =  (2a  +  1)  or  m  =  (2a  +  3). 

The  remainde”  of  the  proof  is  identical  with  the  preceeding,  and  will  be  omitted 
to  avoid  repetition. 

4.  Balanced  Designs  for  the  2m  Factorial  Experiment 

Several  balanced  fractions  are  available  for  the  2m  factorial  experiment 
which  permit  estimates  of  main  effects  free  of  two  factor  interactions.  Certain 
fractions  will  now  be  presented  along  with  a  note  on  their  optimality. 

Def.  4.1.  If  T1  and  T^  are  two  balanced  fractions  with  N  assemblies  each,  the 
efficiency  of  T^  relative  to  T?  will  be  defined  as 

(4.1)  Eff(T  /T  )  =  [tr(EE’):1][tr(EE’)“i]_1 

2 


Th?  efficiency  Eff  T),  of  a  balanced  fraction  T,  is  defined  as  the  relative 
efficiency  of  T  relative  to  a  (possibly  nonexistent)  fraction  satisfying  the 
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optimality  condition  of  theorem  3.4.  We  now  present  some  actual  designs  obtainable 
-by  the  method  of  the  last  section. 

Series  I.  Let  T  =  [J  :  0  :  I  :  I  3,  where  I  is  obtained  from  the  matrix  I 

— -  mp  mp  m  m  m  m 

by  interchanging  0  and  1.  for  these  fractions,  j T j  =  2(m+p),  and 


Eff(T) 


4[m(nHp )  -  4(m-l)][2m  +  p  -  2 j 


(m+o)~  (2m+P-l)  (m+p-1)  [(m+p)(m-l)  -  4(m-2)3 
where  c  equals  0  or  1  according  as  (nup)  is  even  or  odd.  It  can  be  checked  by 
direct  calculation  that  the  efficiency  of  these  designs  decreases  as  m  and  p 
increase,  from  about  70%  for  m  =  7,  o  ^  o,  to  about  25%  for  m  =  12,  p  =  5. 

As  an  example  of  a  fraction  constructed  from  balanced  incomplete  block 
designs,  consider  the  series  of  8IBD  with  parameters  (v*  -  m  =  ^A*  +  3,  b*  = 

4A*  4  3,  y*  =  2X*  +  1,  K*  =  2A*  +  1,  A*),  where  4A*  +  3  is  a  prime  power, or  a 
prime.  Their  existence  is  proved  in  Bose  (1940).  Let  be  the  incid®'  •  matrix 
of  a  BIBD  of  this  form. 

Series  II.  Let  T  =  u’  @  0  ,  ST,  ®  T  .  Then  N  =  4(2A*  +  2)  and  u  =  (2A*  +  2). 

The  number  of  d.f,  for  error  is  1.  Also,  the  efficiency  of  these  fractions 

equals  1,  so  that  they  are  optimal.  It  is  interesting  to  note  that  if  the  two 
assemblies  represented  by  J  ^  and  0  ^  are  omitted  from  the  design,  the  efficiency 
is  reduced  to  4(A*  +  1)/(8A*  +  5). 

Series  III.  Let  2A*  +  1  <  m  <  4A*  +3,  and  cut  out  (4A*  +  3  -  m)  varieties  from 

the  BIBD  given  in  series  II.  Then  a  design  with  parameters  (m,  b  =  4A*  +  3, 

r  =  2A*  t  1,  A)  is  obtained,  and  with  to  =  (2A*  +  1).  From  theorem  3.5  the  optimum 
value  of  w  is  (2A*  +2),  and  the  efficiency  of  the  design  may  be  computed  by 
substitution.  If  a  column  of  i’s  and  a  column  of  0's  is  adjoined  to  T,  then  N  = 

4(2A*  +  2)  and  to  =  (2A*  +  2)  and  the  efficiency  of  the  designs  in  this  series  becomes  1. 
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In  the  following  table  we  give  fractional  designs  of  resolution  IV  for 
the  2m  series  when  7  <_  m  12,  and  for  several  practical  values  of  N.  To  obtain 
the  set  of  assemblies  in  the  design,  first  write  the  incidence  matrix  T*  of  a 
BIBD  with  parameters  as  given  in  the  table,  and  then  cut  out  any  v*  -  m  rows. 

One  such  BI3D  in  each  case  is  given  in  Cochran  and  Cox  [(I960),  pp.  469]  with 
plan  number  as  given  in  the  table.  Then  adjoin  (N/2  -  b)  columns  of  zeros  to  T* 
and  call  this  matrix  T^.  The  final  set  of  assemblies  is  then  T  =  T^  @  T^  where 
TV.  is  obtained  from  T^  by  interchanging  1's  and  0's.  It  should  be  noted  that 
there  are  several  BIBD's,  in  addition  to  the  one  listed,  which  will  give  fractional 
designs  with  the  same  parameters.  From  the  point  of  view  of  efficiency  these 
designs  are  equivalent;  however,  they  may  differ  with  respect  to  the  number  of 
degrees  of  freedom  for  error.  This  aspect  of  the  designs  is  too  cumbersome 
to  be  studied  here,  except  that  we  have  (N  -  2m)/2  <_  ng  N  -  2m  by  Theorem  3.4. 

The  efficiency  given  in  the  table  for  the  designs  is  relative  to  the  class 
of  all  balanced  fractional  designs.  Within  the  class  of  designs  of  the  form 
Tf  ®  T^,  where  T,  is  the  incidence  matrix  of  a  BIBDU,  the  designs  given  are 
optimal. 


(b*r*tf) 

£ 

Plan 

(7,3,1) 

1.00 

7 

(7,3,1) 

1.00 

7 

(7,3,1) 

0.88 

7 

(11,5,2) 

0.94 

19 

(11,5,2) 

1.00 

19 

(13,4,1) 

1.00 

22 

(14,7,3) 

0.93 

10 

(15,6,2) 

0.98 

16 

(16,6,2) 

1.00 

27 

(15,7,3) 

1.00 

25 

(16,6,2) 

0.95 

27 

(19,9,4) 

0.99 

31 

(18,8,3) 

1.00 

11 

(11.5.2) 

(11.3.2) 
(13,4,1) 


(15. 

,6,2) 

(16 

,6,2) 

(16 

,6,2) 

(16 

,6,2) 

(19 

,9,4) 

(19 

,9,4) 

19 
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then  we  can  use  only  16  of  them,  since  16  is  the  multiple  of  8  nearest  co  but  not 
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greater  than  22.  This  would  result  in  a  loss  of  6/22  of  the  available  information. 

Thus  the  purpose  of  this  paper  is  to  obtain  good  nonorthogonal  or  irregular 
designs  of  resolution  IV  for  the  2  series.  Besides  being  of  desirable  size*  a 
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the  estimates.  In  this  paper,  such  designs  (with  even  N)  are  obtained,  'fiiese 
designs  are  optimal  with  respect  to  the  trace,  determinant  and  the  largest  root 
criteria  which  are  shewn  to  be  equivalent.  In  other  words,  among  all  possible 
designs  a  given  value  of  N,  our  designs  minimize  the  trace,  the  determinant  and 
the  largest  root  of  V. 
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